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We present both experimentally and theoretically the transformation of radially and azimuthally 
polarized vector beams when they propagate through a biaxial crystal and are transformed by the 
conical refraction phenomenon. We show that, at the focal plane, the transverse pattern is formed 
by a ring-like light structure with an azimuthal node, being this node found at diametrically opposite 
points of the ring for radial/azimuthal polarizations. We also prove that the state of polarization of 
the transformed beams is conical refraction-like, i.e. that every two diametrically opposite points of 
the light ring are linearly orthogonally polarized. 

OCXS: (260.5430) Polarization; (260.1440) Birefringence; (260.1180) Crystal optics 


I. INTRODUCTION 

Vortex beams with radial and azimuthal polarizations 
are the paradigm of vector beams [1]. The nature of the 
electric field make them ideal candidates for applications 
where extremely tightly focused beams are needed, such 
as in material processing [2] or super-resolution imaging 
[3]. However, the propagation of radially and azimuthally 
vortex beams through anisotropic media such as uniaxial 
or biaxial crystals remains almost unexplored [4-8]. 

In anisotropic crystals, the most particularly spectac¬ 
ular phenomenon is conical refraction (CR), occurring in 
biaxial crystals [9-11]. In CR, a homogeneously circu¬ 
larly polarized focused Gaussian beam passing along one 
of the optic axis of a biaxial crystal is transformed into 
a pair of concentric bright rings with a dark (Poggen- 
dorff) ring, at the focal plane of the system [9-12], see 
Fig. 1(a). The state of polarization (SOP) of the rings is 
the following: any point of the rings is linearly polarized 
with the azimuth 4> rotating along the ring so that every 
two diametrically points are orthogonally polarized, i.e. 
4> = , where ip is the azimuth in cylindrical coordi¬ 

nates and Pc is the orientation of the plane of optic axes 
of the crystal [11, 13]. By moving symmetrically from 
the focal plane in the axial direction, the CR transverse 
pattern becomes more involved and two intensity peaks, 
known as Raman spots [14], are found at axial positions 

2 : given by 2:Raman = ±^J\poZR [15], where zr is the 
Rayleigh range of the focused input beam. Therefore, 
the CR beam forms an optical bottle closed axially by 
the two Raman spots, as shown in Fig. 1(b). Note that 
all these features of the CR beam are found for a circu¬ 
larly polarized input beam of Gaussian transverse profile 
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satisfying the condition pQ = Rq/wq ^ 1, which covers 
most of the applications of CR [15-19]. For a linearly po¬ 
larized input beam, the CR ring possesses an azimuthal 
node corresponding with the polarization of the ring with 
opposite azimuth to the input beam, see Fig. 1(c). 


'-15 -10 -5 0 5 10 15 

z/Zr 

FIG. 1. (Color online) Intensity distribution of a CR beam 
under conditions of po = Rq/wq ^ 1. (a) Bright and dark 
Poggendorff rings at the focal plane obtained with a circularly 
polarized Gaussian input beam. Blue double-arrows indicate 
the polarization distribution along the rings, (b) Intensity 
distribution in the z-r plane showing the free space evolution 
of the CR beam, (c) CR transverse intensity pattern at the 
focal plane for a horizontally polarized Gaussian input beam. 

CR related to radially and azimuthally polarized 
beams, in what follows named as R/A beams, has been 
studied recently showing that R/A beams can be ob¬ 
tained from CR by using an interferometric system [20], 
by using linearly polarized Bessel beams [21] and from 
the transformation of input R/A vortex beams by means 
of CR [22]. In this work, we use a similar experimental 
set-up as in [22] to show the beam evolution and SOP of 
the CR beams obtained beyond the biaxial crystal when 
input R/A vector beams are used (Section II). In Sec¬ 
tion III by means of the diffractive theory of CR [9, 10, 23] 
we perform the numerical simulations corresponding to 
our experimental scheme, finding a reasonable agreement 
between both numerical simulations and experimental re¬ 
sults. Finally, arguments to understand the origin of the 
differences between our results and the ones reported in 










2 


[22] and the summary of our work are provided in Sec¬ 
tion IV. 


II. EXPERIMENTAL PROCEDURE 
A. Experimental set-up 

The experimental set-up is sketched in Fig. 2. A diode 
laser at 640 nm coupled to a monomode fiber generates 
a Gaussian beam that passes through a linear polarizer 
(LP), a quarter waveplate (QWP) to ensure circular po¬ 
larization and a polarization state generator (PSG). The 
PSG consists in 12 linear polarizers joint to form an az- 
imuthally segmented polarizer that transforms the input 
polarization into radial, see top part of Fig. 2. Therefore, 
after passing through the PSG a Gaussian R-polarized 
beam is obtained. To generate an A-polarized beam, two 
half wave-plates (HWPs) with a mutual angle of 45° are 
additionally arranged into the system. Unfortunately, 
the two HWPs used do not operate perfectly for the wave¬ 
length used and we have checked that the A-polarized 
beam obtained has an ellipticity of nearly 3°. The R/A 
beam is focused with a lens of 200 mm focal length and 
passes through a biaxial crystal and parallel to one of its 
optic axis. The biaxial crystal is placed always before the 
expected focal plane of the beam in the absence of the 
crystal. An additional imaging lens is used to transfer 
the transverse planes along the beam propagation direc¬ 
tion onto the GGD camera. A LP is used as analyzer, to 
check the state of polarization of the beam after being 
transformed by the GR phenomenon. We use a commer¬ 
cially available (GROptics) KGd(W 04)2 biaxial crystal 
with a = 16.9 mrad and length I = 23 mm yielding GR 
ring radius of Rq ~ 390 jam. For more details about the 
material and the alignment procedure see e.g. Ref. [11]. 


B. Experimental results 

Fig. 3 presents the experimental transverse intensity 
patterns at the focal plane an at the Raman spots ob¬ 
tained from a R-polarized (first row) and an A-polarized 
(second row) input beam, in the absence of the linear 
analyzer in Fig. 2. At the focal plane and for both in¬ 
put polarizations the transverse pattern is analogous to 
what would be obtained in case of a linearly polarized 
input beam: the intensity of the GR rings is azimuthally 
crescent, with a point of null intensity (in the case of an 
input A-polarized beam, i.e. Fig. 3(e), this null intensity 
becomes a region of minimum intensity due to the non¬ 
perfect operation of the half waveplates) . However, by 
comparing Figs. 3(b) and (e) it can be observed that the 
node for an A-polarized input beam appears at a diamet¬ 
rically opposite position along the ring with respect to a 
R-polarized input beam. Note that the maximum inten¬ 
sity points on the transverse patterns at the focal plane 
are found at the position of the azimuthal sector of the 
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FIG. 2. (Color online) Experimental set-up. A circularly po¬ 
larized Gaussian beam is obtained from a diode laser couple 
to a mono-mode fiber, a linear polarizer (LP) and a quarter 
waveplate (QWP). A radially (R) polarized beam is generated 
when the Gaussian beam passes through a segmented radial 
polarizer. To obtain an azimuthally (A) polarized beam, we 
use additionally two half waveplates with 45° rotation of the 
mutual orientation of their axis. The R/A beam is then fo¬ 
cused by a lens (FL, 200 mm of focal length) being its prop¬ 
agation direction paralel to one of the optic axes of a biaxial 
crystal (BC) of length I = 23 mm. An additional LP is used 
to check the state of polarization of the output beam. Finally, 
an imaging lens (IL) is used to project the transverse pattern 
along the axial direction into the CCD camera. 


R/A-polarized input beam coinciding with the GR po¬ 
larization. In contrast, the minimum intensity points are 
found at positions where the R/A-polarized input beam 
is orthogonal to the GR polarization. 

At the position of the Raman spots and for both R/A- 
polarized beams, the transverse pattern is similar to the 
one expected for a uniformly linearly polarized Gaussian 
input beam, with the difference that now a node is formed 
at the beam center. This is expected, since at that point 
there is a polarization singularity due to nature of the 
R/A-polarized beams. The reported experimental trans¬ 
verse intensity patterns are not perfectly uniform due 
to the discreteness of the PSG used to obtain the R/A- 
polarized input beams, which is made of 12 individual 
polarizers joint together. 

Fig. 4 shows the transverse patterns of the GR beams 
at the focal plane when a linear polarizer (LP) is used as 
analyzer and rotated at angles [0,45,90,135]°. Both for 
R- and A-polarized beams, the node introduced by the 
LP appear at the same azimuthal position. Additionally, 
the LP only introduces one node, being this placed at 
opposite positions along the ring when comparing the re¬ 
sults corresponding to 0° and 90° [Figs. 4(a,e) and (c,g), 
respectively], and 45° and 135° [Figs. 4(b,f) and (d,h), 
respectively]. Therefore, we can conclude that the state 
of polarization of the GR is the same for both R/A input 
polarizations and that it is GR-like. 


III. THEORETICAL MODEL 


In order to corroborate to experimental results pre¬ 
sented above, we have performed the corresponding nu- 
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FIG. 3. Experimental transverse intensity patterns obtained 
from a R/A-polarized input beam (first/second row) at (a,d) 
the first Raman spot plane, (b,e) the focal plane and (c,f) the 
second Raman spot plane. 
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FIG. 4. (Golor online) Experimental transverse intensity pat¬ 
terns obtained from a R/A-polarized input beam (first/second 
row) at the focal plane when a linear polarizer (LP) used as 
analyzer after the crystal is oriented at (a,e) 0°, (b,f) 45°, 
(c,g) 90° and (d,h) 135°. 


merical simulations of transformation of R/A-polarized 
Gaussian beams by means of CR. We consider R/A- 
polarized input beams with transverse electric field 


ER = e ^ (cos 0, sin 0)6*"^, (1) 

Ea = (sin0, — cos(/))e*'^, (2) 

where r = and (f) = tan“^(^/x); and Fourier 

transform A{k) 

— Axi^iAjCx AyiyiAjCy^ (3) 

oo 

Jj E:,{r)e~^'^"^dxdy, (4) 

— OO 
OO 

^?/(«) = Jj Ey{f)e-^^-^dxdy, (5) 

— OO 


where R = Note that the spiral phase term 

(exp(i0)) in Eqs. (1) and (2) comes from the generation 


of the R/A-polarization with the segmented polarizer il¬ 
luminated with circular polarization. Beyond the crystal, 
the electric field of the transformed beam is given by the 
following set of equations [23]: 




•Bo,a(f,/9o) = JJ j^^-^Sm{poK)Aa{K)dB) 

— OO 

,a(^ Po) = JJ (cos (pok) + sin (/Oqk)) Aa{K)d^) 


where a = {x^y}^ n is the mean refractive index of the 
biaxial crystal, r = (x,^), k, = and ^ = R-r — 

‘ The expressions for the displacement vector D in 
terms of Eqs. (6) and (7) are 

Bx — B^^yiv , Po) T B\^x{x 5 Po)? (3) 

Dy = Bo,x{e Po) + Bi^y{r, -po). (9) 

Fig. 5 presents the numerically calculated transverse 
intensity pattern of the CR beam for radial and az¬ 
imuthal input polarizations. The corresponding electric 
field in the x and y direction of at the focal plane is 
shown in Fig. 6. As it can be appreciated, the numerical 
simulations obtained with Eqs. (l)-(9) reproduce the ex¬ 
perimental observations shown above: (i) the transverse 
pattern at the focal plane is formed by an annular light 
structure with an azimuthal node, (ii) the position of 
the node is at diametrically opposite points when R/A- 
polarized beams are used, (iii) the state of polarization 
at the focal plane is CR-like, i.e. every two diametri¬ 
cally opposite points are orthogonally polarized; (iv) far 
from the focal plane a vortex-like beam is observed at the 
beam center and (vi) the beam evolution is non symmet¬ 
ric with respect to the focal plane, at variance with CR 
beams obtained from homogeneously polarized Gaussian 
input beams. 


IV. DISCUSSION AND CONCLUSIONS 

In this section we discuss the similarities and differ¬ 
ences of our results with the work of Grant et al [22]. On 
the one hand, the transverse pattern that we have found 
at the focal plane for R/A-polarizations is the same as 
in Ref. [22]. On the other hand, we have found addi¬ 
tionally a vortex-like structure far from the focal plane, 
at variance with what was reported in Ref. [22]. Note 
that a R/A-polarized beam possesses an optical singu¬ 
larity in the electric field at the beam center and this 
optical singularity must be preserved along its transmis¬ 
sion through a linear system, such as a biaxial crystal, 
as we have found. Additionally and what is more impor¬ 
tant, the state of polarization of the CR pattern at the 
focal plane that we have found is CR-like. In contrast, 
in Ref. [22] CR patterns with R/A-polarizations were re¬ 
ported. As shown in [13, 19, 24, 25], the polarization dis¬ 
tribution at the focal plane for po ^ 1 is always linear at 
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To sum up, we have reported the transformation of ra¬ 
dially and azimuthally polarized beams by a biaxial crys¬ 
tal throughout the conical refraction phenomenon and we 
have obtained a reasonable agreement between experi¬ 
ments and numerical simulations. We have demonstrated 
that the transverse intensity pattern at the focal plane is 
formed by two concentric light rings with an azimuthal 
node. The state of polarization of these rings has been 
shown to be CR-like, i.e. such that every two opposite 
points of the rings are orthogonally polarized. Finally, 
we have found that far enough from the focal plane along 
the axial direction a vortex-like beam is obtained and a 
non-symmetric beam evolution. 


FIG. 5. Numerically calculated transverse intensity patterns 
obtained by using Eqs. (l)-(9) from a R/A-polarized input 
beam (first/second row) at (a,d) the first Raman spot plane, 
(b,e) the focal plane and (c,f) the second Raman spot plane. 
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EIG. 6. Numerically calculated transverse intensity patterns 
of the (a,c) Ey and (b,d) Ex fields obtained by using Eqs. (1)- 
(9) from a R/A-polarized input beam (a,b/c,d) at the focal 
plane. 

any point of the ring being every pair of diametrically op¬ 
posite points orthogonally polarized with independence 
of the state of polarization of the input beam. There¬ 
fore, even for non-homogeneously R/A-polarized beam, 
the polarization distribution of the CR rings at the focal 
plane must be CR-like. 

Finally, note that in our case the input beam is R/A- 
polarized with a spiral phase term due to the helicity of 
the circularly polarized Gaussian beam passing through 
the segmented polarizer. In contrast, in Ref. [22] the in¬ 
put beam passing through the biaxial crystal is a perfect 
R/A-polarized beam, since it is generated with a S-plate. 
However, we have numerically checked that all the differ¬ 
ences commented above with respect to the features of 
the CR beams reported in Ref. [22] are still found, see 
Appendix A. 


Appendix A: Numerical calculations for radially and 
azimuthally polarized beams obtained from a S-plate 

In this Appendix we present the transverse intensity 
patterns for the electric field behind the biaxial crystal 
obtained when radially and azimuthally (R/A) polarized 
beams generated with a S-plate are considered. Such 
beams can be described as follows: 

E'ji = e~^ (cos 0, sin 0), (Al) 

E'y^ = e~^ (sin0, — COS0). (A2) 

Fig. 7 presents the numerically calculated transverse in¬ 
tensity pattern of the CR beam for R/A input polariza¬ 
tions. The corresponding electric field in the x and y 
direction of at the focal plane is shown in Fig. 8. As it 
can be appreciated, the polarization distribution of the 
CR beam at the focal plane is as reported in Fig. 6, i. e. 
CR-like. Additionally, a node at the Raman spot can be 
also found, as in the case of R/A-polarized beams with a 
spiral phase term, see Fig. 5. 
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